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Abstract 

Kinetic equations for the process of multiphoton bremsstrahlung of electron in a matter have been used to derive 
the equations for the moments of distributions of electrons over the number of emitted photons and over the energy 
loss. The equations for the moments have been solved in continuous slowing down approximation. It has been 
shown that the moments of multiphoton distributions can be expressed in terms of the moments of macroscopic 
differential cross section of the process. The results of analytical calculations have been compared with the data of 
Monte Carlo simulation. 



1 Introduction 

The bremsstrahlung of electrons passing through a matter is an example of stochastic process where the number 
of emitted photons as well as their energies are random. The importance of multiphoton effects was underlined in |Q] 
where the problems of experimental study of bremsstrahlung cross section were discussed. The summary energy of 
emitted photons, as a rule, is measured in such experiments and the distribution of electrons over energy loss differs 
from the energy distribution of radiated photons Q . Kinetic equations for the distribution of electrons over the number 
of emitted photons and for the distribution over the radiation energy loss are given in the paper. They are used to 
get the equations for the moments of the distributions. Qualitative analysis of the path and energy dependences of 
the moments is given in continuous slowing down approximation with simple relativistic Bethe-Heitler formula for the 
macroscopic differential cross section of bremsstrahlung. It is shown that the moments of differential cross section can 
be expressed in terms of the moments of energy loss distribution which can be determined in experimental way. The 
results of the analysis are compared with the data of Monte Carlo simulation. 



2 Kinetic equations for multiphoton bremsstrahlung 

Let us consider an electron of energy E traveling in a matter the path length I. The probability P(n\E, I) to emit 
n photons of bremsstrahlung along the pass I obeys the balance equation (the Kolmogorov-Chapman equation) |J: 

P{n\E,l) = (l-sE(E))P(n\E,l-s) + sE(E) [ S ^ 7 j,f ^ P(n - l\E - E 1 , l)dE 1 , (1) 

Jo 

T,(E) and Y,(E 1 ;E) being the total and differential macroscopic cross sections of bremsstrahlung, s is a small part 
of I . The first term in the right side of Eq.(^J) corresponds to electrons which pass small path s without radiation 
and 1 — sYj{E) is corresponding probability. These electrons have to emit n photons along the rest of the path I — s. 
The second term corresponds to the electrons which emit the first photon passing path s and s'S(E) is corresponding 
probability. These electrons have to emit n — 1 photons after that but the energy of electrons equals to E — £" 7 where 
E 1 is random energy of the first emitted photon and £(-E 7 ; E)/Y>(E) is the probability density function of E 7 . In the 
limit s — > the formula (Q) gives the integro-differential equation for P(n\E, I) Q|: 

8 f E 

—P{n\E,l) + H{E)P(n\E,l)- £(£ 7 ; E)P(n - l\E - E 1 , l)dE 1 (2) 

with boundary condition 

P(n\E,l = 0)=S no . 

In similar way one can obtain the balance equation for the probability density function P(Q\E,l) describing the 
distribution of electrons over the energy loss Q : 
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P(Q\E, I) = (1 - a Z(E))P(Q\E, l-a) + aE(E) J ^^- P(Q - E 7 \E - E 7 ,l)dE^, 
which leads to the equation j|] 

—P(Q\E,l) + X(E)P(Q\E,l)- / Z(E^,E)P(Q- E n \E - E 1 ,l)dE 1 (3) 

Jo 

with boundary condition 

P(Q\E,l = Q) = 6(Q), 

S(Q) being the Dirac S -function. 

Multiplication of both side of Eq.(||) by Q" and integration over Q gives the equations for the moments of the 
distribution P(Q\E,l) : 

W(E,l)= [ Q n P(Q\EJ)dQ. 
Jo 

The equations for two first moments have the form 

-Q{E,l) + H(E)Q(E,l)- J X(E J ;E)Q(E-E y ,l)dE 7 = (3(E), (4) 

Q(E, I = 0) = 0, 



di 



di 



Q*(E,l) + Z(E)Q2{E,l) - X(EY,EjQ2(E-E 1 ,l)dE 1 ='y(E)+2 E 1 Y I (E 1 -E)Q(E-E 1 ,l)dE 1 , (5) 

Jo Jo 



Q 2 (E,l = 0) = 0. 

Quantities (3(E) and "f(E) in (f|) and (g) are the moments of differential cross section H(Q;E) 



(3(E) = / E 1 Y 1 (E 1 -E)dE 1 , 
Jo 

j(E) = / E 2 H(E 1 -E)dE 1 . 
Jo 



3 Continuous slowing down approximation 

If differential cross section E(2? 7 ; E) is rapidly decreasing function of the variable E 1 the collision integrals in (j|), 
(IH), (j|) can be transformed by the Taylor expansion of integrands: 



P(n - 1\E - E^l) « P(n - I) - E- ( -^P(n - 1\E, I), 



Q(E - B 7 , a Q(£, - E^Q(E, I), 



f(E - E y ,l) « Q2 (i?j ,) _ Ei JL Q 2 {Ej iy 



Such transform leads to the equations: 



-P(n|S, + E(£)P(n|£, J) = £(£)P(n - 1\E, I) - (3(E)— P(n - 1\E, I), (6) 
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^Q(E,l)+l3(E)^Q(E,l)= 0(E), (7) 

j^{E, I) + p{E)-^W{Ei I) = i{E) + 2(3(E)Q(E, I) - 2 7 (£)^Q(£, 0- (8) 
It follows from the Eq.(|6|) that in S(-E) ~ const approximation P(n\E, I) is the Poisson distribution: 

P(n\E, I) = exp(-lE(E)) ^ E ^ U (g) 

TV. 

with the mean number (multiplicity) of emitted photons equals to 

n(E, I) = Y,™ =0 nP(n\E, I) = IZ(E). (10) 
To get the solutions of (0), (||) quantities 0(E) and "f(E) have to be determined. 

4 Numerical data for bremsstrahlung 

For a qualitative analysis of the multiphoton bremsstrahlung we used the simplified relativistic Bethe-Heitler 
formula |§] for the differential cross section: 

nE 1 -,E) = ^(\(l-^) + § U E^<E 1 <E), (11) 

(c) 

where t is the radiation length of the material passed by electron and Ey — aE is the cut off energy. The cut of for 
the bremsstrahlung photon spectra is due to the density effect H which suppresses significantly the soft part of the 
spectrum. Parameter a can be estimated by the formula 

hojp 
a ~ 



where Uuj p is the plasmon energy of the target material, mc 2 — 511 keV is the rest energy of an electron. 
Such Yj(E 1 ;E) leads to the energy independent cross section 

T,(E)^-(\\n--\) (12) 
t 6 a d 

and to the approximate formulas for (3 and 7 : 

m « f , 

F 2 

One can see from (|To|) , jl^ ) that the photon multiplicity n doesn't depend on the electron energy and is determined 
by the target thickness only. 

Substitution of approximate expression for (5 and 7 into (0) , (||) and solution of these equations gives the formulas 
for the moments: 

Q(E,l)=E(l-exp(- 1 -)), (13) 
W(E, I) = ^(3 - 4exp(-|) + cM-j)) 

and for the variance: 

a 2 Q (E,l) = Q 2 (E 7 l)~Q 2 (EJ) = —(4cxp(--) -3exp(-|) - 1). (14) 

It is seen from (|T^), (Q), that for ~ -C 1 the moments of bremsstrahlung differential cross section (3 and 7 can be 
expressed in terms of the moments of the energy loss distribution Q and Oq : 
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Q(E,l)*lj = l{3(E), 

which can be determined in experiments where the characteristics of energy loss distribution are measured. 

Formulas (|Io|), ( |l3|) and (jlj) were used to calculate W, Q and <Jq for 8 GeV electron bremsstrahlung in tungsten 
target. The results were compared with the data of Monte Carlo simulation using the GEANT code J(|. The Monte 
Carlo method was used for calculation of distribution of electrons over the number of emitted photons. The calculations 
show that the distribution strictly follows to the Poisson law @. Mean number of photons (Fig.l) agrees with (|lo|), dl2| ) 
if a ~ 1.1-10 -5 . This value of a corresponds to the cut of energy E)f* » 88keV. Fig. 2 shows that the energy distribution 
P(E 7 ;E, I) of emitted photons (the number of photons per unit interval of E y ) really vanishes for E 7 < E^ c \ 




Fig. 1: Mean number of photons. Fig. 2: Energy distributions of emitted photons at 

solid line — Eq. (10), points — simulation. 1=0. 01cm 

For small I when multiphoton effect can be neglected the energy distribution P(E 7 ; E, I) can be expressed in terms 
of the differential cross section H(E^;E) : P(E 7 ; E,l) — IY,(E 7 ;E). Fig. 3 shows that the distribution P(E 7 ; E,l) 
obtained by simulation agrees with that one calculated with using the formula (p"l|). 

Data on Q(E,l) are given in Fig. 4 as well as the mean energy of emitted photons: E 7 = 2. It is seen that E 7 
slowly decreases with I because of electron energy decreasing. The variance of energy loss (Tq(E, I) is shown in Fig. 5. 
One can see that the analytical formulas ( |l3|) , ((TJ) for Q and Oq agrees with the data of Monte Carlo simulation quite 
good. 

The Monte Carlo code was used to calculate the probability density function P(Q\E 7 1) describing the distribution 
of electrons over energy loss Q and the ratio of P(Q\E, I) to the energy spectra. The results are given in Fig. 6 
which displays the difference between radiation loss spectrum from finite thickness target and intensity spectrum of 
bremsstrahlung in the field of single nucleus. It is clear that multiphoton emission depresses the low energy part of 
the distribution and increases the high energy part. The effect rises with the target thickness increasing. 

This work is supported by the Grant of the Ministry of General and Professional Education of the Russian Feder- 
ation. 
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Fig. 3: Comparison of photons energy distribution cor- 
responding to Eq. (11) (solid line) with simulation at 
I = 0.01cm 



Fig. 4: Mean energy loss Q(E, I) and 
mean energy of emitted photons E 1 . 
Solid line — Eq. (13), points — simulation 
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Fig. 5: Variance of energy loss. Solid line — Eq. (14), points — simulation 
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Fig. 6: The ratio V (Q, I) = P g°Jg|2, 
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